Spinon confinement and the Haldane gap in SU(n) spin chains 
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We use extensive DMRG calculations to show that a classification of SU(«) spin chains with regard to the 
existence of spinon confinement and hence a Haldane gap obtained previously for valence bond solid models 
applies to SU(«) Heisenberg chains as well. In particular, we observe spinon confinement due to a next-nearest 
neighbor interaction in the SU(4) representation 10 chain. 
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Introduction. — The properties of quantum spin chains have 
been a vital area of research in condensed matter physics. 
Starting with Bethe's solution of the spin \ Heisenberg model 
(HM) in 1931—, the field quickly evolved and significantly 
influenced many other areas of physics. In accordance 
to previous findings by Andrei and Lowenstein 2 , Faddeev 
and Takhtajan- observed in 1981 by consideration of Bethe 
Ansatz solutions that the elementary excitations of spin chains 
are spinons carrying spin 1/2. From the identification of the 
0(3) nonlinear sigma model as the low-energy field theory 
of antiferromagnetic SU(2) spin chains, Haldane conjectured 
in 1983 that chains with integer spin possess a gap in the 
magnetic excitation spectrum, while a topological term ren- 
ders half-integer spin chains gapless— . This, and in particular 
the gap in the magnetic energy spectrum for integer chains, 
was confirmed by experiment—. An elegant paradigm of the 
gapped spin 1 chain in terms of a valence bond solid (VBS) 
model was given by Affleck, Kennedy, Lieb, and Tasaki 
(AKLT)2. 

Haldane's classification^ applies only to SU(2) spin sys- 
tems. In the field of ultracold atoms and optical lattices, how- 
ever, experimental realizations of SU(«) spin systems may be 
possible in the near future. In particular, it has been proposed 
very recently that SU(«) antiferromagnets with n up to 10 may 
be realized with ultracold alkaline-earth atoms without any 
need to fine-tune any of the interaction parameters^. This per- 
spective makes a general classification of the magnetic spectra 
of antiferromagnetic SU(«) spin chains highly desirable. 

It is, however, not immediately clear how to address the 
question which SU(«) spin chains generally possess a Hal- 
dane type gap, as Haldane's original work cannot directly be 
generalized to SU(n); previous attempts were limited to spe- 
cial cases^i£. One possible route is to interpret the Haldane 
gap as the zero point energy of the oscillator describing the 
relative motion of pairwise confined spinons — , and to em- 
ploy simple paradigms to determine whether spinons in a spin 
chain with spins transforming according to a given represen- 
tation are confined or not. Following this line of reasoning, 
two of us— recently formulated a variety of VBS models 
for SU(n) spin chains, and investigated which models exhibit 
spinon confinement and hence a Haldane gap. This led to a 
classification of SU(«) spin chains into three categories, as re- 
viewed below. The conjecture which motivated this work was 
that the "gap/no gap" classification found for the VBS mod- 
els would apply to general SU(«) spin chain models, and in 
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FIG. 1: (Color online) (a) Spinon confinement in the AKLT model. 
The potential is proportional to the distance \x\ between the spinons 
and allows an interpretation of the Haldane gap as the zero point 
energy of a linear oscillator, (b) Spinons in the spin 1/2 VBS state 
are domain walls between the two ground states and hence free. 



particular to the SU(n) HMs. 

In this paper, we explicitly confirm this conjecture using ex- 
tensive Density Matrix Renormalization Group (DMRG) cal- 
culations. In particular, we observe confinement of spinons 
and hence a change of the universality class triggered by a 
next-nearest neighbor interaction in the simplest example of 
the third category, the SU(4) rep. 10 model. 

Spinon confinement in VBS states. — Arovas et ai— pointed 
out that the spin 1 AKLT model //aklt = Zi^CM + 1) 1S 
a good starting point to reach the Heisenberg point Hum = 
2/5/5 /+ i = 3 }Zi(P2(i, i+ 1) + l/3Pi(/, i+ 1) - 2/3) perturba- 
tively. Here Pg(i, . . . , i+m) is a projector of the total spin of 
m + 1 neighboring sites onto representation R of SU(2). Both 
models show an excitation gap according to Haldane's con- 
jecture. For the AKLT model, exact calculations — showed 
that the static spin-spin correlations decay exponentially, be- 
ing the signature of an excitation gap (see also the second line 
in Tab.[T]i. The AKLT model contains a characteristic length 
scale associated with the correlation decay, which may only 
change quantitatively, but not qualitatively upon moving to 
the HM. In the following, we assume that also beyond SU(2), 
the existence of a length scale at a VBS point persists as one 
moves to the Heisenberg point. This links the gap property of 
a VBS model with the HM for the same spin representation. 

It is easy to illustrate how the Haldane gap can be traced 
to spinon confinement in the AKLT state. It can be written 
as a product of local spin singlets consisting of two "virtual" 
spins 1/2 placed on adjacent sites. Projecting onto the sym- 
metric subspace on each site gives the spin 1 representation: 
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Here, each circle denotes a "virtual" spin 1/2 and each line a 



singlet bond. The dashed box illustrates exemplarily a lattice 
site. The AKLT Hamiltonian is just the sum over the projec- 
tors onto the subspace with total spin 2 of adjacent sites. For 
two neighboring sites of Yaklt, a s P m singlet and two individ- 
ual "virtual" spins 1 /2 are present, i.e., ® i ® | = ffi 1. The 
tensor decomposition of two spin 1 representations, however, 
additionally contains a spin 2 subspace, i.e. , 1 ® 1 = ffi 1 ©2. 
Hence, it follows that ^aklt is annihilated by //aklt, whereas 
all other states will be lifted to finite energy. To create a pair 
of spinons, one breaks a singlet bond (see Fig.QJ). The region 
between the two spinons yields an energy penalty for //aklt, 
or stated differently, the two spinons feel a confining potential 
V(x) = F\x\. As the Hamiltonian of the relative motion of the 
two spinons consists of a kinetic energy and this linear con- 
fining potential, we obtain a linear oscillator. We interpret the 
zero point energy of this oscillator as the Haldane gap. 

For half-integral spin chains, the corresponding VBS is 
the Majumdar-Ghosh (MG) model— (see also the first line 
in Tab.[T]i. In our terminology both spontaneously n-merized 
(like MG) and translational invariant states (like AKLT) are 
VBS states. Its ground state consists of local singlet bonds 
placed on adjacent sites. In contrast to the AKLT model, the 
MG ground state is two-fold degenerate and invariant only 
under translations by two lattice sites. As in the AKLT model, 
elementary excitations are created by breaking a singlet bond. 
The two emerging spinons can serve as domain walls inter- 
polating between the two ground states (see Fig.QJi), i.e., the 
spinons do not feel any confining potential. Accordingly, the 
spin-spin correlations decay abruptly, i.e., <S,S, + 2) = 0, and 
no associated length scale describing an exponential decay 
emerges. Moving from the MG point to the isotropic HM, one 
passes a phase transition from the dimer phase into a critical 
spin liquid phase 15 . The deconfined spinons of the MG model, 
however, remain deconfined when moving to the generic HM 
and give rise to the observable two-spinon continuum—. Ap- 
plied to arbitrary integer and half integer representations, this 
line of argument reproduces Haldane's original conjecture for 
SU(2) spin chains. This interpretation is consistent with ear- 
lier attempts to describe the spin 1 elementary excitation of 
related models as confined solitons— or two-spinon compos- 
ite particles—. 

Previous conjecture. — So far, we reproduced Haldane's 
conjecture for SU(2) in terms of spinon confinement. We 
now review the general classification of the magnetic exci- 
tation spectra for SU(n) chains motivated by the VBS states 
introduced previously— . This classification depends only on 
the degree n of SU(«) and the number of boxes A in the Young 
tableau (YT) associated with the respective spin representa- 
tion (see also Tab. [Q: (I) If A and n have no common divisor, 
the model will support free spinon excitations and not exhibit 
a Haldane gap. (II) If A is divisible by n, the model will ex- 
hibit spinon confinement and hence a Haldane gap. (Ill) If A 
and « have a common divisor different from «, the presence 
of the confining potential will depend on the range of interac- 
tions. If q is the largest common divisor of A and «, interac- 
tions ranging to the (n/q)-th neighbor are required for spinon 
confinement. 

For SU(2), the classification is identical to Haldane's, as the 
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TABLE I: Classification for spin-spin correlations of certain SU(«) 
VBS states as introduced in Ref.|T2|. Listed are the spin representa- 
tion, the corresponding Young Tableau (YT), the spin-spin correla- 
tion length <fi2, and the category C w.r.t. our conjecture. The SU(4) 
rep. 10 VBS is the simplest model of category III. The correlations 
for the MG and the rep. 6 model decay abruptly and no correlation 
length can be defined. 



third category becomes accessible only for n > 4. If this clas- 
sification applies to the HMs as well, as we advocate here, we 
will find gapless spectra, unique ground states, and algebraic 
correlations for the first category. For the second category, we 
will find gapped spectra and exponentially decaying correla- 
tions, as it is possible to write the HM as VBS model plus a 
small perturbation (similar as for the spin 1 case—). For the 
SU(4) rep. 10 chain belonging to the third category, the gap- 
less spectrum of a nearest-neighbor HM should acquire a gap 
as a next-nearest neighbor interaction term is added. 

In Tab. Q] we further list the exact values for the static 
spin-spin correlation lengt h 19 ' 20 ! 21 . We expect that Heisenberg 
chains for all representations without a well-defined VBS cor- 
relation length £ (like line 1 and 3 in Tab. will be gapless. 
For representations with a well-defined correlation length 
we have to distinguish a "standard" gapped category II and the 
non-trivial category III, which we will treat separately below. 

Verification of the conjecture. — To verify our predictions 
for the SU(n) HMs, we performed extensive DMRG 22 - 23 i 24 i 25 
studies for all SU(«) representations up to the ten-dimensional 
representations of SU(3) and SU(4). Note that we use the 
"Abelian" quantum numbers of SU(«), as the "non-Abelian 
version" of DMRG— is not really convenient for SU(3) and 
SU(4) representations. In order to treat the large block sizes 
of our DMRG calculation we use a Posix thread parallelized 
DMRG code 27 . In our list of representations, we find gapless 
models such as the fundamental representations of SU(3) {A = 
1 in YTbox notation), SU(4) (A = 1), and the six-dimensional 
representation of SU(3) (A = 2). The ten-dimensional repre- 
sentation of SU(3) ( i i i i . A = 3) provides an example for a 
gapped model. As the local basis dimension can become very 
large, we are not always able to extrapolate to the thermody- 
namic limit and read off the gap from there. Instead, we use 
the entanglement entropy S a = -Tr[p a logp,,], where p a is 
the reduced density matrix in which all the degrees of free- 
dom on sites i > a are traced out. This way, there is no need 
to extrapolate. In particular, for critical spin models associ- 
ated with conformal field theories (CFT), the entanglement 
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FIG. 2: (Color online) Entanglement entropy (48 sites, PBCs) for the 
critical SU(n) Heisenberg chains with SU(3) rep. 3 (Nan = 5000 kept 
DMRG states), rep. 6 (2V cut = 5000), rep. 10 (AU = 1650), and SU(4) 
rep. 4 (N cut = 12000). Except of SU(3) rep. 10, all models are gapless 
and nicely follow the sinusoidal shape according to Q~|l. The latter 
clearly develops a plateau structure as expected for gapped models. 



entropy is given by— 




where L is the chain length, c\ a non-universal constant, 
and c the central charge of the associated CFT. In case 
of Heisenberg-like SU(«) spin chains, the CFT is a Wess- 
Zumino-Witten model with topological coupling k (SU(«)jt 
WZW) with central charge c = k(n 2 - l)/(k + n). We use 
periodic boundary conditions (PBCs) rather than hard wall 
boundary conditions in our DMRG calculations, as the latter 
show even/odd oscillation triggered by the boundaries. These 
would render the interpretation of our results for the category 
III models below ambivalent. 

If the models exhibit a gap in the excitation spectrum, the 
entanglement entropy does not have the sinusoidal shape but 
saturates after a few sites. If the models have dimerized (or 
trimerized, tetramerized, etc.) ground states, this manifests 
itself in the entanglement entropy by oscillations of the en- 
tropy with an oscillation period of two sites (or three, four, 
etc.) reflecting that the ordered states get pinned by numeri- 
cal noise. The numerical results are shown in Fig. [2] For the 
nearest neighbor SU(«) HMs with fundamental representation 
(A — 1), we find the sinusoidal shape (Q]) with fitting parame- 
ter c - n - 1 for n = 2, 3,4—. For the SU(3) representation 6 
HM, we find a critical model with a central charge larger than 
c = 2, with the difference due to marginal operator perturba- 
tions. For larger chain lengths, the associated inverse logarith- 
mic corrections disappear and the central charge approaches 
c = 2, as expected. The situation is hence comparable to the 
spin 3/2 HM— . This implies that the SU(3) representation 6 
HM belongs to the SU(3)i WZW universality class. For the 
ten-dimensional representation of SU(3) (A - 3), we find that 
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FIG. 3: Entanglement entropy of the SU(4) rep. 10 HM with next- 
nearest neighbor interactions for (a) J 2 = 0.57] and (b) J 2 = 0.1571, 
as discussed in the text. The inset in (a) illustrates the SU(4) rep. 10 
VBS state where each circle denotes a fundamental SU(4) represen- 
tation and the horizontal line denotes antisymmetric coupling. Four 
connected circles hence represent an SU(4) singlet. 



the entanglement entropy saturates after a few sites. This in- 
dicates a gap (see Fig.|2|. 

Category III models. — We have further performed extensive 
calculations for the HM of the ten-dimensional representation 
of SU(4) (A = 2) with nearest and next-nearest neighbor in- 
teractions. This is the simplest example belonging to the third 
category. Let us first look at the corresponding VBS state 
shown in the inset of Fig. [3^. The state is two-fold degen- 
erate and invariant under translations by two lattice spacings. 
Its parent Hamiltonian involves three-site interactions^. The 
correlations decay exponentially with length £ (see Table|I|. 

For nearest neighbor interactions, we find that the model is 
gapless and critical as required by the Affleck-Lieb theorem^. 
The fitted central charge is larger than 3 in accordance with 
the SU(4)i WZW universality class in the presence of loga- 
rithmic corrections. Following our categorization, applying a 
next-nearest neighbor interaction should force the system in 
a gapped phase with confined spinons. As a first test, we ap- 
proximate the Hamiltonian of the VBS state within the J\—Jj 
model //vbs ~ Si^A+i + l/2S;Si+2- The corresponding en- 
tanglement entropy is shown in Fig. [3^. It shows a strongly 
oscillating behavior indicating dimerization, as expected from 
the analogy with the VBS state shown in the inset. The con- 
jecture that the HMs behave as the associated VBS models 
thus is confirmed for this type of models. It still remains to 
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TABLE II: Critical couplings 15 i 30 i 31 i 32 of J { -J 2 models of different 
SU(w) reps, are shown. For representations belonging to the cat- 
egory II with Haldane gap, no phase transition occurs. All other 
SU(ra) models tend to rc-merize (i.e., dimerize for n = 2, etc.)—&. 
We further conjecture, that the critical couplings {JilJi)c increase 
when n and/or A is increased. Consequently, the SU(4) rep. 10 chain 
is expected to tetramerize for a critical coupling which is larger than 
(J 2 /Ji)c ~ 0.5. For the SU(3) rep.6 and the SU(4) rep.4 model, only 
preliminary results are available, as indicated by the asterisks. 

be shown, however, that this gap is not due to a tetramer- 
ized phase as expected for a frustrated SU(4) model driven 
by Ji- To rule this out, we have collected the critical cou- 
plings of all relevant SU(«) /1-/2 models in Tab.HIl From 
there, one would expect that a transition into a tetramerized 
phase for SU(4) rep. 10 would only occur for some critical 
coupling (J2/Ji) c > 0.5. By contrast, we find a dimer phase 
rather than a tetramer phase, and for couplings smaller than 
the expected critical coupling. In our DMRG calculations, we 
observe the existence of the dimer phase for a value as small 
as J2 ~ 0. 15/i (see Fig.[3ji). The oscillation amplitude of the 
dimerization becomes weaker but is still present. Our data 



cannot be understood in the context of frustrated HMs, but 
are fully consistent with our conjecture of spinon confinement 
triggered by a next-nearest neighbor interaction. 

Physically, the opening of a gap can be understood as fol- 
lows. As we apply a next-nearest neighbor coupling J2, pairs 
of neighboring sites effectively cluster into new sites, and 
pairs of spins transforming under the original representation 
form new representations on the new sites. The state hence 
dimerizes on the original lattice, while it remains invariant 
under translations of the new lattice. The original represen- 
tations implied that spinons were deconfined and no gap oc- 
curred. Spinons in the effective model with spins transforming 
under the new representations, however, are confined and the 
spectrum is gapped. The next-nearest neighbor coupling J2 in 
our category III model here hence produces an effect similar 
to the effect produced by a coupling of two spin 1/2 chains 
into a spin ladder. 

Conclusion. — We have shown that a classification of SU(«) 
spin chains with regard to the existence of a Haldane gap ob- 
tained previously for VBS models applies to SU(«) Heisen- 
berg chains as well. The results provide evidence in favor 
of our hypothesis that this gap can be interpreted as the zero 
point energy of an oscillator describing the relative motion of 
confined spinons. 
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